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ABSTRACT 


The  techniques  for  producing  n+-p-n+  structures  in  gallium 
arsenide  are  described.  Evaluation  of  the  surface  state  densities  of 
anodic  oxide  on  gallium  arsenide  show  minima  at  about  mid  gap  of 
2 x 10u  eV  cm2.  These  densities  are  affected  by  annealing  conditions 
and  their  measurement  by  sweep  time. 

Recommendations  for  future  work  include  the  use  of  indium 
phosphide  as  the  semiconductor. 


***** 


1.  INTRODUCTION 

(a)  The  objective  of  the  program  is  to  develop  vertical  chan- 
nel insulated-gate  field  effect  transistor  devices  in  GaAs  which  will  be 
ultimately  capable  of  delivering  5 watts  of  power  (Class  A)  with  minimum 
power  gain  of  6 dB  over  the  frequency  range  4-8  GHz.  Additional  design 
goals  for  these  devices  are  that  they  shall  be  linear  in  phase  (±5°  devi- 
ation) and  gain  (IMD  <-20  dB)  over  the  range  of  operation. 

(b)  Review  of  Previous  Accomplishments 

The  previous  accomplishments  of  this  program  have  been  recorded 
in  detail  elsewhere, ^ but  it  is  appropriate  to  review  the  main  areas  of 
progress  in  the  fabrication  of  a vertical  MOS  GaAs  FET. 

1 . 1 Device  Fabrication 

Devices  have  been  fabricated  using  the  planar  geometries  shown 
in  Figure  1.  The  devices  have  shown  FET  action  as  illustrated  in 
Figure  2.  The  source  and  drain  contacts  for  these  devices  is  gold  ger- 
manium with  an  anodic  oxide  insulator  (1000X  thick)  and  an  aluminum  gate. 
The  ohmic  source  and  drain  contacts  are  alloyed  for  10 secs  at  4S0°C  in 
argon  containing  10%  hydrogen. 


H 
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1 . 2 Etch  Development 


As  part  of  the  program,  etches  were  developed  to  produce  the 
V-groove  structure  (shown  in  Figure  3)  as  part  of  the  processing 
required  for  the  V-groove  device  which  will  be  discussed  later.  They 
consisted  of  mixtures  of  the  form  X : H2O2  : H2O  where  X can  be  HC1, 

HF  or  NaOH.  The  ratio  (in  mis)  of  the  components  x : y : 100  where  x, 
the  amount  of  X (from  2 to  6 mis)  does  not  have  much  effect  on  the  etch 
rate  and  y,  the  amount  of  H2O2  determines  the  etch  rate  from  the  rela- 
tionship 

Etch  rate  = 483. 3y  A/min  . 

This  relationship  holds  from  2 to  8 mis  of  H2O2.  At  10  mis  and  above 
the  surface  of  the  gallium  arsenide  begins  to  develop  a matt  finish. 

The  HF-based  etch  is  anisotropic  and  causes  rounded  mesas, 
whereas  the  NaOH  etch  is  isotropic,  resulting  in  the  structures  shown 
in  Figure  4. 

Anodization  etching  has  also  been  used  to  remove  material  and, 
in  addition,  is  able  to  delineate  the  p-n  junction  structures.  This  is 
illustrated  in  Figure  5. 

1.3  Anodization 

The  most  favorable  dielectric  in  this  program  has  been  anodic 
oxide  formed  on  gallium  arsenide  using  a 1 : 1 mixture  of  3%  citric  acid 
and  ethyleneglycol  buffered  with  ammonium  hydroxide  to  a pH  of  5.5.  This 
solution  had  the  minimum  self-etching  rate  (0.024  A/sec ^)  and  pro- 
duced uniform  films  with  good  dielectric  properties  (breakdown  field  = 
4.37  x 106  V/cm,  e /e  = 7.2).  The  anodization  procedure  is  compatible 
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with  device  fabrication  and  the  anodizing  solution  does  not  attack  any 
of  the  device  structures. 

1.4  Proposed  Device  Structure 

(a)  V-groove  Structure 

The  V-groove  structure  is  planned  as  a stepping  stone  between 
the  planar  devices  and  the  final  vertical  channel  VMIST  geometry 
(described  below).  The  reasons  for  using  this  geometry,  shown  in 
Figure  6,  is  the  greater  accessibility  for  inspection  than  the  VMIST 
device  and  the  somewhat  simpler  fabrication  procedure.  This  fabrication 
procedure  is  outlined  in  Table  1.  The  mask  set  for  these  devices  is 
shown  in  Figure  7 and  the  fabrication  procedure  in  Figure  8. 

Initially,  the  structure  for  the  V-groove  devices  will  be  an 
n«--p-n+  structure  as  shown  in  Figure  6(a).  Ultimately,  there  will  be 
an  additional  n-  layer  as  shown  in  Figure  6(b)  for  the  realization  of 
the  VMIST  geometry. 

An  analysis  of  the  effect  of  photolithographic  design  rules  on 
the  power  density  of  the  V-groove  device  is  given  in  Appendix  A. 

(b)  VMIST  Structure 

The  VMIST  structure  is  the  form  of  the  device  ultimately 
intended  for  this  program  and  follows  closely  that  already  developed  in 
silicon. ^ A cross-section  of  the  device  is  shown  in  Figure  9. 
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Fig.  6-layer  structures  for  V -groove  (a)  and  VMIST  (b)  devices 
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Figure  7 Mask  set  for  V-groove  devices 
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TABLE  1 

Fabrication  of  V-groove  Device 

1)  Evaluate  n+-p-n+  starting  material  for  layer  thicknesses. 

2)  Sputter  Si(>2  (0.5  urn  thick)  onto  surface  wafer. 

3)  Photoresist  - Mask  #724-1  to  delineate  device  mesas. 

4)  Using  HF  : H2O2  : H2O  in  the  proportions  1 : 3 (30\)  : 100  (mils) 
etch  mesas  [Fig.  8(a)]. 

5)  Remove  Si02. 

6)  Sputter  Si02  (0.5  um  thick). 

7)  Open  V-groove  opening  - Mask  #724-5  - with  the  required  opening 
depending  on  the  epi  layer  thicknesses. 

8)  Etch  V-groove  using  NaOH  : H2O2  : H2O  in  the  proportions 
1 gm  : 3 mis  (30%)  : 100  mis  [see  Fig.  8(b)]. 

9)  Open  gate  opening  - Mask  #724-4  [see  Fig.  8(c)],  etch  oxide  over- 
hang, anodize  GaAs,  evaporate  gold  gate  metal  and  reject  excess 
with  photoresist. 

10)  Open  source  drain  openings  - Mask  #724-2  - etch  oxide,  remove  resist. 


11)  Photoresist,  delineate  source-drain  contact  openings  - Mask  #724-3  - 
evaporate  gold-germanium  and  reject  excess  metal  [Fig.  8(d)]. 

- 13  - 


1 . 5 Three  Key  Technologies 

(a)  n*-p-n4  Structure 

Of  the  various  approaches  that  are  possible  to  fabricate  the 
structure,  the  following  were  chosen  for  reasons  of  availability  and 
control . 

(i)  Liquid  Phase  Epitaxy.  Liquid  phase  epitaxial  layers 

2 an  thick  of  p-type  (N^  = 8 x 1016  cm'3)  material  have  been  grown  on 
n*  substrates  (N^  = 1017  cm"3).  On  top  of  this,  a further  n-layer  with 
Np  * 10 17  cm'3  and  thickness  0.5  ym  has  been  grown. 

(ii)  Ion  Implantation.  Ion  implantation  of  silicon  in  VPE 
p-type  layers  on  n-type  substrates.  The  carrier  concentrations  of  the 
VPE  p-type  layers  are  shown  later  in  Table  2.  The  thickness  of  the 
p-layer  shown  in  Figure  6(a)  can  be  relaxed  to  a value  as  high  as  5 microns 
during  initial  ion  implantation  experiments  to  allow  for  evaluation  of 

the  layer  and  some  freedom  in  the  implant  energies  and  annealing  condi- 
tions. Excess  material  can  be  readily  removed  by  chemical  etching  or 
anodizing  etching. 

(b)  Stable  MDS  Oxide 

fl  2) 

Previous  workv  ’ J has  indicated  that  anodic  oxides  held  great 
promise  for  GaAs  MOS  devices.  Anodization  is  a very  convenient  process 
to  incorporate  into  a device  fabrication  procedure  since  it  is  very  con- 
trollable and  is  compatible  with  other  fabrication  steps.  Our  results 
indicate  that  hysteresis  is  still  a problem  as  far  as  long-term  stability 
of  the  device  is  concerned. 
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( 2 ) 

Results  obtained  on  previous  planar  devices  have  shown  FET 


action  but  there  is  no  positive  evidence  that  inversion  of  the  p-type 
surface  occurred.  Referring  to  Figure  10,  it  is  likely  that  n-type 
material  was  left  between  the  n-mesa  on  the  p-type  substrate  and  that 
accumulation  of  electrons  occurred  in  the  depleted  n-region. 


- IS  - 
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Anodic  Oxide 


Thin  Depleted 
N-Type  Layer 


Fig.  10-Cross  section  of  planar  device  showing  possible 
presence  of  thin  n -layer  under  gate 
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2.  MATERIALS  DEVELOPMENT 


2.1  LPE  Techniques 

Liquid  phase  epitaxy  (LPE)  was  employed  to  produce  five  layers 
suitable  for  the  vertical  transistor.  One  of  the  major  strengths  of  LPE 
is  the  ability  to  quite  precisely  control  both  n-  and  p-type  doping. 
Furthermore,  since  the  p-type  layer  was  constrained  to  be  5 x 1016  cm-3, 
this  makes  VPE  a less  attractive  choice  than  LPE  vis-a-vis  p-doping 
control. 

Figure  11  shows  a picture  of  the  LPE  apparatus  and  Figure  12 
shows  the  graphite  boat  used  to  grow  multilayer  structures.  In  order  to 
achieve  an  LPE  abrupt  junction  vertical  structure  with  good  doping  con- 
trol, it  is  necessary  to  choose  dopants  that  have  small  segregation 
coefficients  and  are  slow  diffusers.  For  the  n-type  layer  Sn  was  chosen 
and  Ge  was  selected  for  the  p-type  dopant.  Both  are  Group  IV  elements 
and  are  amphoteric;  thus,  the  effective  segregation  coefficient  (i.e., 
the  n-  or  p-typeness)  is  a strong  function  of  temperature.  Figure  13 
shows  the  segregation  behavior  of  Sn  as  a function  of  temperature.  For 
our  particular  device  a growth  temperature  of  800°C  was  selected  to  aid 
in  substrate-melt  wetting  which  gives  good  surface  morphology  for  device 
fabrication.  Experiments  were  undertaken  to  determine  the  segregation 
coefficient  of  Ge  at  800°C.  The  segregation  coefficient,  Kge,  is  defined 
as  follows: 


(na-nD) 


R.T. 


Ge 


Atomic  Concentration  of  Ge  in  Melt 
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Figure  11  Liquid  phase  epitaxial  growth  system 


RM-71450 


Tr|5H 

(•  — ‘ ■ 

luff 

Grams  of  Sn 


Curve  684984- C 


(N0  -N^l  cm 


-3 


Figure  13  Grams  of  Sn  needed  to  achieve  particular  doping  level  in 
5 gm  Ga  melt  with  growth  temperature  as  a parameter 


- 20  - 


Our  data  yields  a value  for  K_  of  3.49  x 10" 3 for  LPE  growth  in  the 


vicinity  of  800°C. 

Growth  conditions  were  set  for  initial  melt  saturation  of 
800°C.  After  a melt  bake  of  approximately  1 hour,  a 12#C/hr  cooling 
ramp  was  initiated.  After  1°C  of  growth  on  a dummy  substrate,  the  seed 
crystal  was  slid  into  place  to  grow  the  p- layer  for  6 min.  At  this 
point,  the  slider  was  adjusted  so  that  the  seed  was  removed  from  contact 
with  the  melt  and  the  dummy  substrate  was  placed  under  the  n-melt. 

After  S min  of  growth  from  the  n-melt  on  the  dummy  substrate,  the  seed 
was  positioned  under  the  n-type  melt  and  growth  proceeded  for  3 min. 
Using  this  technique,  we  were  able  to  provide  layers  with  the  following 
characteristics: 


p-type  2p  8 7.7  x 1016 

n-type  lw  8 1.5  x 1017 

In  addition,  thicx  p-type  layers  with  7.7  x 1016  net  hole  concentration  1.4 

were  supplied  for  ion  implantation  experiments. 

ifl 

2.2  Vapor  Phase  Epitaxy 

Vapor  epitaxial  growth  techniques  have  been  developed  to  prepare 
all  epitaxial  multilayer  GaAs  VMIST  structures  and  epitaxial  ion  implanted 
structures  according  to  the  following  schemes: 


all  epitaxial  epi  ion  implanted 

VMIST  Structures 


ion  implanted 
p 5 um 
n+  substrate 


n+  1 pm  1018  cm-3 
p 1 pm  mid  1016  cm"3 
n“  3-4  pm  < 1015  cm"3 
n+  substrate 


Epitaxial  layers  were  deposited  on  polished  silicon  doped  n+ 
substrate  wafers  orientated  with  100  surfaces.  Substrate  wafers  were 
obtained  from  Laser  Diodes,  Inc.  and  subsequently  polished  in  an  acid- 
peroxide  solution.  Epitaxial  layers  were  grown  using  the  AsCl3/Ga/H2 
technique  within  the  multilayer  epitaxial  deposition  system  shown  in 
Figure  14. 

The  n~  drift  region  was  grown  undoped  with  net  carrier  concen- 
tration in  the  1013  cm'3  range  using  high  AsCl 3 mole  fraction  growth 
conditions.  An  example  of  a 1.8  yin  thick  undoped  epitaxial  layer  con- 
centration profile  on  an  n+  substrate  is  shown  in  Figure  15. 

P-type  epitaxial  layers  were  grown  zinc-doped  by  addition  of 
a dilute  vapor  of  diethyl  zinc  in  hydrogen  to  the  reactor  growth  stream. 
The  diethyl  zinc  was  obtained  in  5 nines  purity  from  Alpha  Products  and 
thermostated  at  -20°C  to  give  suitable  vapor  pressure  for  doping.  Molar 
flow  rates  of  diethyl  zinc  calculated  to  be  in  the  10'6  moles  min-1 
range  gave  p-type  epitaxial  depositions  of  6 x 1015  to  8 x 1017  cm-3. 

The  actual  concentration  depended  on  growth  rate  and  mole  fraction  of 
AsCl 3 used  in  the  deposition.  The  optical  micrograph  of  Figure  16  shows 
a cross-section  of  a 1.7  ym  thick  p surface  layer  and  a 5 ym  thick  n~ 
drift  layer  on  an  n+  substrate. 

Sulfur  doped  n+  surface  layers  of  about  1018  cm'3  were  grown 
by  dilute  H2S  in  H2  to  the  growth  stream  and  from  a heated  solid  sulfur 
source.  Attempts  were  made  to  grow  n+  layers  using  a liquid  tin  chloride 
source,  but  the  high  vapor  pressure  over  tin  chloride  proved  difficult  to 
control,  and  this  technique  was  abandoned.  A multilayer  (n+)  - (p)  - (n-)  - 
(n+  substrate)  epitaxial  structure  is  shown  in  Figure  17  with  a 2.8  ym 
n+,  a 3.5  ym  p,  and  a 2.8  ym  n'  drift  region. 
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Concentration,  cm 


Single  p-type  epitaxial  layers  were  grown  by  diethyl  zinc 
doping  for  ion  implantation  and  are  typified  by  the  6.3  \m  thick  layer 
shown  in  Figure  18.  Concentrations  of  p-type  epitaxial  layers  were 
evaluated  by  zero  bias  capacitance  measurements  on  aluminum  Schottky 
barrier  contacts. 

A compilation  of  zinc-doped  single  and  multilayer  epitaxial 
growth  wafers  is  shown  in  Table  2.  The  semi- insulating  substrate  sam- 
ples were  used  for  evaluation  purposes. 


6.3  ym 


CD4-R25-2  1700X 


Fig.  18  - 6.3  um  p layer  grown  in  n-type  substrate. 


RM- 76841 


TABLE  2 

Zn-doped  Single  and  Multilayer  Epitaxial  Growth  Wafers 


Wafer  No. 

Substrate 

n" 

D 

n+ 

p Concentration 

CD3-R40-1 

n+ 

4.2  urn 

-2 

n+ 

2.1 

6 x 1015  cm-3 

CD4-R32-1 

n+ 

7.0 

9 x 1016 

-2 

n+ 

2.5 

6 x 1016 

CD4-R33-1 

n+ 

2.6-  3.5 

1.3  x 1017 

-2 

n+ 

- 

- 

CD4-R34-1 

n+ 

5.2 

6 x 1017 

-2 

n+ 

3.1 

8 x 1017 

CD4-R35-1 

n+ 

9.1 

3.7  x 1017 

-2 

n* 

4.2 

4.9  x 1017 

CD4-R36-1 

n+ 

2.8 

-2 

n* 

2.8 

- 

CD3-R38-1 

n+ 

2.8 

3.5 

-2 

n+ 

1.5 

1.7 

CD3-R39-1 

n+ 

3.8 

0.5 

-2 

n* 

1.4 

<0.5 

CD3-R41-1 

SI 

5.2 

3.5  x 1016 

-2 

SI 

2.1 

1.4  x 1016 

CD3-R42-1 

SI 

4.2 

8.4  x 1017 

-2 

SI 

8.5  x 1017 

CD3-R43-1 

SI 

5.0 

4.2 

-2 

SI 

6.0 

2.0 

CD4-R22-1 

SI 

12.6 

6.4  x 1017 

-2 

SI 

9.6 

4.4  x 1017 

CD4-R23-1 

SI 

2.0 

6 x 1016 

-2 

SI 

0.5 

7 x 1016 

CD4-R24-1 

SI 

5.5 

-2 

SI 

6.0 

CD4-R25-1 

SI 

6.3 

6 x 1016 

-2 

SI 

2.8 

5 x 1016 

CD4-R28-1 

SI 

3.1 

-2 

SI 

1.4 

CD3-R36-1 

SI 

8.4 

0.7 

5 x 1017 

-2 

SI 

5.0 

0.5 
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2.3  Ion  Implantation 

Silicon  ions  were  implanted  into  several  of  the  epitaxial 
layers  shown  above.  The  slice  numbers  and  the  required  concentrations 
and  thicknesses  are  shown  in  Table  3.  The  implants  were  made  with  a 
full  dose  of  Si*  at  150  keV  and  a quarter  dose  at  50  keV. 


TABLE  3 

Samples  to  be  Implanted 


Slice  No. 

REQUIRED 

Concentration 

Layer  Thickness 

CD4-R32-2-1  1 

CD4-R33-1-A 

CD4-R35-1-1 

6 x 1017 

0.2  um 

CD4-R3S-2-A 

CD4-R32-2-2 

CD4-R33-1-B 

> 1.2  x 1018 

0.2  urn 

CD4-R35-1-2 

CD4-R35-2-B  , 

A calculated  profile,  according  to  LSS  theory,  is  shown  in 
Figure  19  for  the  first  group  of  slices  shown  in  Table  3.  The  lower 
energy  provides  a flat  carrier  profile  to  the  surface. 

Measurements  using  a four-point  probe  were  made  on  Sample 
CD4-R33-1-A  implanted  with  1 x 1013  at  150  kV  and  2.5  x 1012  at  50  kV. 
With  no  annealing  the  sheet  conductivity  o^  ranged  between  3.53  to 


i 

i 


i 

. 
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5.02  x 10“6  ft'1  cm'1  in  the  dark  at  1.5  volts  applied  and  had  a value 
of  3.1  x 10~5  R'1  cm”1  in  the  dark  with  22.5  volts  applied.  This 
implies  that  there  may  be  some  contact  problems  at  the  lower  voltage. 

With  the  application  of  a bright  tungsten  light,  the  lower  voltage  con- 
ductivity values  increased  to  5.7  x 10” 6 ft"1  cm'1;  while  for  the 
22.5  volt  supply  there  was  evidence  of  breakdown  to  the  lower  n-type 
layer  (o^  * 6.7  x 10” 5 O'1  cm'1). 

The  sample  was  then  cleaved  in  two,  and  the  two  halves  desig- 
nated a and  b.  CD4-R33-1-Aa  was  capped  with  600X  of  pyrolytic  S^Ni, 
(deposited  at  680°C)  covered  by  2000X  of  pyrolytic  SiC>2  (deposited  at 
420°C) . A low  temperature  of  750°C  was  chosen  for  the  annealing  to  avoid 
diffusion  of  the  zinc  from  the  p-type  layer  as  much  as  possible.  The 
annealing  was  performed  in  forming  gas.  The  values  of  sheet  conductivity 
after  annealing  are 

3.67  x 10'**  ft'1  cm'1  (at  1.5  volts  in  the  dark) 

5.17  x 10'**  ft'1  cm'1  (at  22.5  volts  in  the  dark) 

5.45  x 10~4  ft'1  cm'1  (at  22.5  volts  in  tungsten  light) 

These  values  should  be  compared  with  a sheet  conductivity  of  6 x 10' 3 ft'1  cm'1 
which  would  be  obtained  if  the  implanted  silicon  were  100%  activated  and 
the  resulting  carriers  had  a mobility  of  3000  cm2/V-sec.  This  figure  for 
mobility  is  pessimistic  (for  1.25  x 1013  dose  at  400  kV  other  workers^ 
have  obtained  figures  closer  to  4000  cm2/V-sec).  Obviously,  the  activa- 
tion is  not  100%,  and  if  we  assume  the  above  mobility,  the  activation  is 
quite  low  - more  like  9%. 


The  second  half  of  Sample  CD4-R33-1-A  was  capped  with  pyro- 
lytically  deposited  Si02  on  the  front  only.  A 650°C  anneal  for  1/2  hour 
yielded 

= 2 x 10*4  (at  22.5  volts  in  dark) 

Ojj  = 3-4  x 10'4  (at  22.5  volts  in  light) 

which  is  somewhat  worse  than  the  standard  encapsulant.  650°C  was  chosen 
both  for  the  zinc  diffusion  problem  and  the  fact  that  gallium  is  known 
to  diffuse  through  Si02. 

Both  samples,  when  examined  using  the  thermal  probes,  showed 
p-type  behavior  (a)  being  more  strongly  p-type  than  (b) . While  this 
measurement  is  not  conclusive,  it  confirms  the  poor  activation  efficiency 
of  the  silicon  implant  into  the  p-type  epitaxial  layer.  Capacitance 
voltage  measurements  (shown  in  Figure  20)  also  indicate  that  the  surface, 
at  least,  is  still  p-type. 

1 
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3.  DIELECTRIC  INVESTIGATIONS 


3. 1 Discussion 

The  two  most  important  factors  that  govern  the  performance  of 
inversion  mode  field  effect  transistors  are  (a)  the  surface  mobility  and 
velocity  of  the  carriers  in  the  inversion  layer  (in  this  case  electrons). 
If  the  mobility  it  not  sufficiently  higher  than  that  of  silicon 
(200  cm^/volt  sec)  then  the  advantage  of  gallium  arsenide  will  have 
disappeared  and  the  more  difficult  technology  development  cannot  be 
justified.  A factor  of  ten  (2000  cm^/volt  sec)  is  a reasonable  target 
for  which  to  aim.  The  saturated  surface  velocity  should  exceed  that  of 
silicon  (5  x 10^  cm/sec)  by  a factor  of  two  (see  Appendix  A)  if  improved 
device  performance  is  expected,  (b)  The  electrical  and  chemical  stability 
of  the  dielectric  used.  Chemical  stability  requires  that  the  layer  not 
deteriorate  on  contact  with  the  air,  be  unduly  hygroscopic  and  be  able 
to  withstand  any  subsequent  fabrication  steps,  annealing  or  bonding 
operations.  Electrical  stability  is  more  subtle  than  the  chemical 
stability.  It  requires  that  the  dielectric  strength  of  the  layers  not 
deteriorate  but,  more  importantly,  that  any  surface  states  present  be 
minimal  and  unchanging. 

Experience  with  silicon  indicates  that  there  are  two  kinds  of 
states  which  affect  MOS  performance:  (a)  Fast  surface  states  (Ngj)  which 
have  communication  with  the  conduction  band  of  the  semiconductor  and 
which  are  filled  and  emptied  as  the  Fermi  level  at  the  surface  is  raised 


and  lowered.  The  effect  of  these  states  is  to  distort  the  shape  of  the 
C-V  curve  as  shown  in  curve  B of  Fig.  21.  Curve  A is  the  ideal  curve 
without  surface  states; 


(b)  Insulator  charges  (Qss)  which  are  trapped  in  the  insulator 
itself.  An  example  of  such  insulator  charges  are  sodium  ions  in  silicon 
dioxide  on  silicon  devices.  These  are  generally  known  as  slow  states 
because  of  their  long  time  constants.  Their  effect  is  to  shift  the 
curve  from  B to  C. 

3.2  Stability 

(a)  Analysis  of  the  Results 

The  major  problem  in  evaluating  the  C-V  data  obtained  from 
measurements  at  1 MHz  of  the  anodic  oxide  or  gallium  arsenide  is  the 
presence  of  the  hysteresis.  A typical  curve  obtained  fc - p-type  gallium 
arsenide  is  shown  in  Fig.  22.  Beginning  at  a positive  voltage  (usually 
15-20  volts)  at  point  A,  the  bias  is  decreased  until  it  reaches  a 
negative  voltage  at  point  B (VA  = 1VB) . During  this  excursion,  the 
semiconductor  goes  from  "inversion"  at  point  A to  accumulation  at  point  B 
following  the  curve  through  point  C.  The  voltage  is  held  at  point  B for 
a fixed  time  and  then  the  voltage  is  returned  at  the  same  rate  to  point  A 
but  this  time  the  curve  proceeds  via  point  D. 

The  theoretical  curve  for  the  same  carrier  concentration  and 
oxide  thickness  but  with  no  surface  or  oxide  charges  present  is  shown  as 
the  dotted  line  in  Fig.  22.  The  measured  curves  are  obviously  shifted 
toward  the  negative  voltages  corresponding  to  the  trapping  of  positive 
charge  but  in  addition  to  the  shift,  the  shapes  of  the  curves  are  perturbed 
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Fig.  22-C-V  curves  showing  hysteresis  loop 


from  the  theoretical  shapes  indicating  the  presence  of  fast  states  whose 
occupancy  changes  during  the  bias  sweep. 

The  hysteresis  is  due  to  trapping  centers  in  the  forbidden 
band  filling  up  as  the  voltage  is  driven  from  A to  B and  then  held  at 
B.  In  order  to  investigate  these  centers,  attention  was  centered  on 
the  curve  ACB  from  whose  shape  the  distribution  of  the  states  in  the 
bandgap  and  their  density  can  be  determined. 

If  the  shift  of  the  capacitance  voltage  curve  from  the  ideal, 
trap  free  case  is  AV  then  the  amount  of  charge  Q trapped  in  the  system 
that  produces  this  shift  is  given  by  Q = CQ  • AV  where  CQ  is  the  capacitance 
of  the  anodic  oxide  layer.  The  surface  potential  of  the  oxide  (<J>S)  can  be 
determined  from  the  measured  capacitance  of  the  system  (C)  so  that  by 
measuring  AV  as  a function  of  C,  we  can  calculate  the  total  amount  of 
trapped  charge  as  a function  of  surface  potential.  By  differentiating 
the  charge  as  a function  of  surface  potential,  we  can  determine  the 
density  of  surface  states  as  a function  of  their  position  in  the  energy 
bands. 

The  derivation  of  the  equations  and  the  computer  program  based 
on  these  equations  is  given  in  Appendix  B. 

The  curve  from  D to  E shows  the  deep  depletion  effects  in  the 
semiconductor  when  the  voltage  sweep  is  too  fast  to  allow  a sufficient 
number  of  minority  carriers  to  be  thermally  generated  to  maintain 
equilibrium. 
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3.3  Effect  of  Treatment  of  Surface  State  Density 


(a)  Annealing  Time 

O 

Fig.  23  shows  the  C-V  curves  obtained  on  a 1000A  thick 
(50  volt)  anodic  oxide  on  p-type  (10^  cm-^)  gallium  arsenide  annealed 
for  various  times  in  a nitrogen  atmopshere  at  360°C.  Using  the  analysis 
described  above  the  surface  state  density  as  a function  of  surface 
potential  is  shown  in  Fig.  24. 

The  minima  of  the  curves  are  located  at  approximately  -0.3  eV 
(on  the  conduction  band  side  of  the  mid  gap  point  which  is  0.0  eV  on 
this  figure) . The  curves  are  not  smooth  in  these  measurements  due  to 
the  fact  that  the  data  was  taken  from  the  C-V  curves  plotted  on  an  XY 
recorder.  Subsequent  measurements  are  taken  directly. 

The  obvious  conclusion  to  be  drawn  from  this  data  is  that 
the  minimum  value  of  surface  density  decreases  as  the  time  of  annealing 
is  increased.  However,  the  curve  obtained  from  the  four  hour  anneal, 
while  having  a lower  minimum,  lies  above  the  3,  2 and  1 hour  curves 
in  the  region  from  -0.85  to  0.2  eV.  Thus  there  is  an  increase  in  the 
surface  state  density  in  this  region.  This  is  borne  out  by  the  size 
of  the  hysteresis  loop  for  the  four  hour  C-V  curve  shown  in  Fig.  23. 

The  rapid  rise  of  the  surface  state  density  as  the  surface 
potential  departs  from  the  position  of  the  minimum  means  that  it  is  very 
difficult  to  shift  the  Fermi  level  very  far  from  this  value.  The  fact 
that  the  barrier  heights  of  Schottky  contacts  on  gallium  arsenide  are 
not  strongly  dependent  on  the  metal  used^  confirms  this  situation. 
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Curve  697429-A 


Fig.  23-CV  curves  as  a function  of  annealing  time 


Surface  State  Density  I cm  eV 


Curve  697433-A 


Sample  P-13  P-Type 

Na=2.15x  1017  cm-3 
920  A Thick  Anodic  Oxide 
Annealed  in  N2  at  360 °C 
Scan  Speed  0.1  V/ sec  r 


I* 


o No  Anneal 
o 1 Hour  Anneal 

* 2 Hours 
■ 3 Hours 

• 4 Hours 


Surface  Potential  ( eV) 


Fig.  24-Surface  state  density  as  a function  of  surface  potential 
for  various  anneal  times 


(b)  Sweep  Time 


The  effect  of  sweep  time  on  the  measurement  of  the  surface 
state  density  is  shown  in  Figure  25.  In  this  sample  the  minimum  lies 
at  near  the  mid  gap  point  in  contrast  to  the  previous  sample.  This 
shift  in  the  position  of  the  minimum  is  not  understood. 

At  slow  scan  speeds  (0.01  volts/sec)  the  minimum  of  the  curve 
lies  at  5.5  x 1011  eV'1  cm-2,  while  at  1 volt/sec  it  is  down  to  1.1  xlO11. 
This  effect  is  related  to  the  hysteresis  loops  seen  on  all  the  eV  mea- 
surements in  that  the  states  have  a finite  charging  up  time  associated 
with  them. 

The  shape  of  the  two  curves  is  very  similar,  indicating  that 
the  time  constant  of  the  traps  is  the  same  irrespective  of  their  posi- 
tion within  the  forbidden  band. 
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Valence  Band 


Fig.  25-Apparent  surface  state  density  variations  with 
sweep  rate  of  voltage 


4.  DEVICE  DESIGN  CONSIDERATIONS 


4.1  Mask  Set  for  V-groove  Device 

The  mask  set  was  discussed  in  Section  1.4(a).  This  device 
does  not  have  an  optimized  design  for  high  frequency  performance  since 
there  is  significant  overlap  between  the  gate  and  the  source  and  drain 
regions  (see  Figure  8).  The  parasitic  capacitances  between  the  gate  and 
drain  is  a feedback  element  that  reduces  the  device  gain. The  device 
is  useful  in  that  it  does  not  require  critical  control  of  the  etching 
procedure. 

The  performance  of  the  V-groove  device  with  optimum  parameters 
is  discussed  in  more  detail  in  Appendix  A. 


5.  FUTURE  PLANS  IN  PURSUING  VERTICAL  CHANNEL 
MOSFET  DEVICE  IN  GaAs 


5. 1  VMIST  Geometry 

This  geometry  has  already  been  proven  in  silicon^3-*  and  is 
the  preferred  geometry  for  this  contract.  Unfortunately,  problems  with 
material  preparation  have  not  allowed  even  the  simpler  V-groove  device 
to  be  fabricated.  In  any  future  effort  the  control  of  the  layer  forma- 
tion would  have  to  be  sufficient  that  the  device  fabrication  procedure 
described  in  Section  1.4(a)  could  be  implemented. 


5.2  Stable  Insulator 

Anodic  oxides  have  proved  interesting  but  have  not  provided 
sufficient  stability  and  freedom  from  hysteresis  to  make  them  a com- 
pletely satisfactory  gate  insulator.  Recent  results  by  Chang  and 
(8) 

others'-  1 have  indicated  that  oxidation  of  the  gallium  arsenide  surface 
using  plasma  oxidation  and  including  a thin  (lOOA)  layer  of  AI2O3  pro- 
duces a dielectric  with  low  leakage  and  a fixed  (slow)  charge  density 
of  5 x 1010  cm*2.  The  film,  however,  does  have  hysteresis  and  no  infor- 
mation was  given  regarding  the  fast  states. 


5.3  n^-p-n*  Profile 

For  the  fabrication  of  the  VMIST  device,  precise  control  of 
the  thicknesses  of  the  various  layers  is  required.  Hie  optimum  way  of 
producing  the  material  would  seem  to  be  the  vapor  phase  or  liquid  phase 
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epitaxy  of  the  n"  and  p layers  and  rely  on  the  ion  implantation  of 
silicon  for  the  n*  layer  to  form  the  source  contact  of  the  device. 

The  layers  discussed  in  Section  2.2  were  fabricated  using 
the  VPE  reactor  shown  in  Figure  26  which  was  designed  to  provide  flexi- 
bility in  doping  of  the  epitaxial  layers.  The  loping  schemes  used  were 
dilute  H2S,  solid  sulfur,  tin  chloride,  thin  doped  gallium  source  (all 
n-type)  and  liquid  diethyl  zinc  (p-type). 

Inadequate  control  of  gas  flow  and  gallium  source  kinetics 
have  led  to  unreproducible  layers  so  that  a new  VPE  reactor  has  been 
designed  (see  Figure  27). 

In  order  to  accurately  control  the  gas  flow,  mass  flow  con- 
trollers were  installed,  and  to  cut  down  on  operator  error  or  variations, 
switching  of  gases  is  accomplished  with  solenoid  valves.  The  mass  flow 
controllers,  valves  and  thermally  controlled  bubbles  for  ASCI3  and  liquid 
dopants  are  all  contained  within  a nitrogen  enclosure  (see  Figure  28)  to 
prevent  possible  oxygen  contamination. 

To  improve  the  gallium  source  kinetics  the  temperature  of  the 
source  has  to  be  accurately  controlled  and,  more  importantly,  the 
temperature  gradient  across  the  source  must  be  close  to  zero  (to  avoid 
partial  crusting  of  gallium  arsenide  on  the  gallium  source).  This  is 
achieved  by  using  a sodium  heat  pipe  24"  long  which  controls  a region 
18"  long  to  an  accuracy  of  temperature  of  +0.5°C  and  a temperature  gra- 
dient less  than  0.5°C  across  the  18"  length.  The  sodium  heat  pipe  and 
the  quartz  deposition  tube  within  it  are  designed  for  2-inch  diameter 
GaAs  substrates. 
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Figure  27  - New  VPE  reactor. 


Figure  28  - Mass  controllers  valving 
and  doping  equipment  in 
nitrogen  enclosure. 


Dual  source  capability  is  provided  on  this  reactor  so  that  we 
may  switch  from  a doped  source  to  an  undoped  one,  and  vice  versa. 

A six- zone  conventional  furnace  provides  a uniform  and  con- 
trollable temperature  profile  over  the  deposition  region. 

The  reactor  to  date  has  been  used  to  grow  unintentionally 
doped  buffer  layers  for  the  ion  implantation  of  the  active  n-type  layer. 
It  has  shown  the  capability  of  growing  good  surface  finish  layers  with 
less  than  10%  variation  in  layer  thickness.  It  has  not  been  used  to 
grow  intentionally  doped  n or  p layers  to  date. 
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6.  FUTURE  RECOWENDATIONS 


The  results  obtained  in  this  program  and  by  other  workers 
indicate  that  the  problem  of  the  large  density  of  surface  states  on 
gallium  arsenide  has  not  been  conquered.  A recent  approach,  plasma 
oxidation  of  the  surface, ^ has  indicated  that  the  fixed  charge  (flat 
band  voltage  shift)  can  be  made  as  low  as  5 x 1010  cm'2.  Devices  have 
been  fabricated  using  anodic  oxidation,  but  most  of  these  are  used  in 
the  depletion  mode  where  surface  states  are  not  a problem  since  they 
aTe  largely  empty.  It  is  not  clear  that  efficient  inversion  of  the 

gallium  arsenide  surface  has  ever  been  achieved. 

(9) 

In  contrast  to  this,  recent  results  have  indicated  that 
indium  phosphide  has  been  inverted  at  the  surface  using  Si02  as  the 
dielectric  layer.  Indium  phosphide  has  a higher  saturated  drift  velocity 
than  that  of  gallium  arsenide  which  leads  to  a higher  frequency  response 
for  a device  with  the  same  geometry.  In  addition,  the  breakdown  field 
of  indium  phosphide  is  higher  than  that  of  gallium  arsenide,  resulting 
in  a higher  operating  voltage  and  hence  higher  power.  These  arguments 
have  been  made  previously for  Schottky  barrier  gate  devices. 

There  is  also  evidence^ ^ that  ion  implantation  of  n+  dopants 
into  indium  phosphide  results  in  higher  activation  than  similar  ones 
into  gallium  arsenide.  Fabrication  of  the  device  should  be  no  more  dif- 
ficult than  for  gallium  arsenide.  Ohmic  contacts  of  gold-germanium  have 
already  been  developed  and  the  problems  of  the  Schottky  barrier  that 
plague  the  MESFET  device  will  not  be  present  in  the  InP  MISFET. 
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APPENDIX  A 


ANALYSIS  OF  THE  EFFECT  OF  LITHOGRAPHIC  DESIGN  RULES 
ON  POWER  DENSITY  FOR  VERTICAL  MIS  GaAs  TRANSISTORS 
HAVING  ONE  MICRON  p AND  n~  REGIONS 


Assuming  that  we  will  fabricate  a vertical  MIS  GaAs  transistor 
having  a p-region  of  one  micron  thickness  and  an  n”  region  of  approxi- 
mately one  micron  thickness  (corresponding  to  device  having  an  output 
transit  time-limited  frequency  of  about  f = *v  16  GHz),  it  is  of 
interest  to  investigate  the  role  of  the  design  rule  of  "a"  microns  such 
that  for  a given  power  level  the  device  power  density  is  not  excessive. 

The  design  geometry  of  the  VMIST  device  is  illustrated  in 
Figure  A-l  where  it  is  assumed  that  two  V-grooves  3p  deep  are  fabricated 
a distance  M apart.  As  can  be  seen  by  the  figure,  we  assume  that  a 

groove  3y  deep  will  cut  through  1)  the  lp  n"  region  with  its  negligibly 

. * 

thin  ion  implanted  n contact,  2)  the  lp  p-region  and  3)  about  lp  into 
the  n++  source  substrate  contact.  Assuming  an  etch  angle  of  60°  (really 
57.7®  for  the  anisotropic  etches  used),  we  find  that 


M -v 


+ 6a 


'v  ^3p  + 6ajp 


The  figure  illustrates  a VMIST  device  with  the  source  as  top  contact  in 
contrast  to  Fig.  8,  but  this  does  not  affect  the  validity  of  the  treatment. 


D vp.  6446A84 


N + Contact 


1 Micron 


1 Micron 


N Substrate 


Fig.  A-l  -Layout  of  unit  cell  of  VMIST  chip 


Calculation  of  Chip  Side  "b"  in  Microns  to  Obtain  a Given  ac  Power 
Output  Pac 

We  assume  that  all  carriers  are  moving  with  an  effective 
saturation  velocity  v After  Siggv  ' 


®m  eox  vsat^ox  ’ 


(A-l) 


where  W is  the  VMIST  oxide  thickness, 
ox 


Let  us  attempt  to  calculate  the  maximum  current  per  unit 
length  I Y by  assuming  the  existence  of  same  maximum  effective  gate 
voltage  V given  by 

Boa 


v = 8 Er_  V 
max  bDq^  ox 


(A-  2) 


where  E* 


= oxide  breakdown  field  (volts/ cm) 

= reduction  factor  which  accounts  for  lack  of  oxide 
integrity,  edge  effects,  surface  state  reduction 
of  effective  field  effect  voltage,  etc. 


Letting  I = g V , we  have 
max  6m  max 


e v 

ox  sat  „„ 

W 8 EBDraax  Wo 

ox 


°r  Jmax  = eox  vsat  EBD  8 * 


(A- 3) 


If  we  take  the  breakdown  field  for  the  gate  insulator  as  (conservatively) 

EfiDmax  ^ 2 x 106  volts/cm;  if  we  assume  8-vO.S  and  we  take  the  effective 

saturation  velocity  anticipated  for  the  GaAs  VMIST  as  v ^107  while  we 

sat 

take  the  effective  saturation  velocity  for  the  silicon  VMOST  as 


v x 106  cm/sec^1^;  and  if,  finally,  we  assume  that  both  the  VMIST 

and  VMOST  will  have  a breakdown  voltage  capability  of  VgD,v20  volts, 
then  we  can  generate  the  following  table  using  Equations  (A-l,  A-2,  A-3) 
and  the  fact  that  the  maximum  quiescent  dc  power  dissipation  for  Class  A 


operation,  Pdcmax»  is  given  by 


Pdcmax  ^ 


^ax  VBD 


(A-4) 


eoxL 


Also  included  is  the  gate  capacitance  per  unit  length  C _ = — ^ , 

gate  WQX 

where  L is  the  channel  width,  assumed  equal  to  ly  as  in  Figure  A-l. 

We  also  calculate  the  intrinsic  cut-off  frequency  of  each  device 

ft  3 *m/2lTCgate* 

In  Table  A-l  the  values  predicted  for  the  VMOST  device 
(g^  = 16  ymhos/y,  Imax  = 160  yA/y)  are  in  excellent  agreement  with 
results  obtained  during  the  VMOST  Contract  No.  N00014-74-C-0012  with 
ONR,  where  Westinghouse  typically  obtained  18  ymhos/y  and  180  mA/y  for 


these  two  quantities. 


This  gives  partial  credibility  to  the  values 


of  gm  = 32  ymhos/y,  Imax  = 330  mA/y,  Pdcmax  *1*5  mW/y  and  ft 'v  16  GHz 
predicted  for  the  GaAs  VMIST,  if  a surface  saturation  velocity  of 
'vlO7  cm/sec  can  be  obtained. 

If  we  assume  a chip  size  of  Bx  By2,  then  we  can  construct  a 
VMIST  device  having  N stripes  of  cross-section  illustrated  in  Figure  A-l. 
Thus,  N would  be  given  by  B/M  or  N * B/(3  + 6a)  from  Figure  A-l.  Each 
stripe  is  B microns  long  and  dissipates  a quiescent  dc  power  of  1.5  mW/y 
from  Table  A-l.  Thus, 

Pdc  * {2  x 1.5  mW/y  x b|  x #/stripes 

- 2 x 1.5  mW/y  x B2  .. 
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■ 


(A- 5) 


COMPARISON  OF  MAXIMUM  CAPABILITIES  OF  VMIST  AND  VMOSTv  ' DEVICES 


If  we  assume  that  the  output  ac  power  P is  50%  of  the  dc 

uC 

quiescent  dissipation,  then  for  a given  required  P (watts) 

3C 


P a 0.5  (dc  quiescent  power) 

»C 


fl.S  x 1Q-3B2 
3 x 6 a 


(A-6) 


Solving  for  the  square  chip  size  necessary  to  achieve  a 
desired  output  power,  P , given  an  obtainable  photolithographic  design 
rule  of  "a"  microns: 


Chip  Side  B = 25  / (3  + 6a) P microns  (A-7) 

SC 

This  equation  is  plotted  in  Figure  A-2. 

The  dc  quiescent  device  power  density,  PD,  is  given  by  the 
division  of  Equation  (A-5)  by  B2  the  device  area,  or 


Power  Density  PD(watts/cm2) 


3 x 105 
(3  x 6a) 


(A-8) 


This  equation  is  plotted  in  Figure  A-3,  where  it  can  be  seen  that  power 
densities  approaching  that  of  IMPATT  diodes  ( -v 105  W/cm2)  can  be  achieved 
in  the  GaAs  VMIST  device  with  lv  design  rules. 


Discussion 

From  Figure  A-2  it  can  be  seen  (point  A)  that  it  would  seem 
possible  to  design  a VMIST  chip  approximately  .008"  x .008"  square 
(excluding  rads,  etc.)  having  the  capability  of  4 watts  rf  at  approxi- 
mately 5-8  GHz  (estimated)  using  only  2p  design  rules  — a rather  loose 
constraint. 


H 
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Curve  697434-A 


APPENDIX  B 


DEVIATION  OF  TUP.  EQUATIONS  FOR  THE  C-V  CURVE 


The  various  forms  of  the  C-V  curve  in  the  absence  of  any 
surface  states  are  shown  in  Figure  B-l.  Curve  ABE  is  the  so-called 
low  frequency  capacitance  curve  where  the  capacitance  measurement  fre- 
quency is  low  enough  that  the  minority  carriers  in  Inversion  can  respond 
to  the  signal.  This  frequency  is  usually  a fow  tens  of  Hz  and  below. 

At  highor  frequencies  the  capacitance  follows  the  ABC  curve,  and  if 
deep  depletion  occurs  the  ABD  curve. 

The  oxide  measurements  aro  taken  by  illuminating  the  sample  at 
point  C to  ensure  equilibrium,  and  then  sweeping  the  voltage  so  that  the 
capacitance  follows  the  curve  CBA. 

To  derive  the  equations  governing  the  shape  of  the  curve  ABC, 
the  low  frequency  equations  are  used  from  point  A to  B and  then  a high 
frequency  approximation  from  B to  C. 

Tho  low  frequency  capacitance  Cj „ is  related  to  the  oxide 
capacitance  Co  and  the  semiconductor  capacitance  via  tho  relationship 


‘LF 


/Co  " 1 + C/t . • 

A'  O 


(B-l) 


whore 


K c 
o o 


ox 


K E 

o 0 


the  diolectric  constant  of  the  oxide  of  thickness  W 


ox 
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The  semiconductor  capacitance 

((sign  <>f  «s) 


Cs  " Si 


sinh|ug  - Upj  ♦ sinh|upj 

F(Us,uf) 


(B-2) 


where  u is  the  surface  potential  and  u_  is  the  Fermi  level  referenced 

' F fa. 


to  the  center  of  the  bandgap  and  is  given  by  Up  = £n 

Ksc0 

Debye  length  capacitance  given  by  Cn  = — = — , and 

u ld 


n^ 


, Cp  is  the 


/K  e kT 

D V IHn-  "5"  is  the  Debye  length 


“I"  i q 

and  n^  is  the  intrinsic  carrier  concentration. 

F(us’Uf)  = |2[cosh(us-uF)  ' cosh(uF)  + ussinh  (UF)] 
The  voltage  on  the  metal  contact  is  given  by 


1/2 


. (B-3) 


' T [US  * (si*n  of  us)  f(v“f] 


CB-4) 


These  equations  are  appropriate  from  A to  B on  the  capacitance- 

voltage  curve.  From  B to  C the  'high  frequency'  capacitance  Cu_  equations 

HF 

apply. 


HF  _ 1 

Co  ' 1 * VCs 


where  C 


K e 
s o 


s W, 


(B-S) 


HF 


and 


r , ) 

"f  , 

X 1 1/2 

u_ 

- Ld  exp  [-Up/2j 

l8up  - 12  + 4 tn  (up  - 2 

5jJ  - 2 exp 

UF  ' T 

. 4 

, 2 ♦ /2/e 

• 

(B-6) 
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An  additional  equation  used  in  the  calculations  is  that  governing  the 
intrinsic  carrier  concentration  n^ 

nA  » / NgNy  exp^-Eg/klj  where  N£  and  Ny 

are  the  densities  of  states  in  the  conduction  and  valence  bands, 
respectively.  It  appears  in  the  program  in  the  form 

n4  = 4.05  x 10 14  x T3/2  x exp|-Eg/2k*rj . (B-7) 

which  is  appropriate  for  gallium  arsenide. 

Iterative  Calculation  of  Carrier  Concentration 

In  order  to  calculate  the  carrier  concentration  of  the  semi- 
conductor from  the  minimum  capacitance,  it  is  necessary  to  use  an 
iterative  process  based  on  the  Taylor  series  since  Equation  (6-B)  gives 
the  capacitance  as  a function  of  Fermi  level  position. 


Let  C ® ^[uf] 

CB-8) 

Define  a function  F(Uf)  * c “ f|up 

(B-9) 

Now  the  Taylor  series  expansion  of  F^upj  is 

F(uf)  ■ F(uF„)*(uF-“Fn)F'  (»F„j 

| ♦ higher  terms. 

(B-10) 

where  upn  is  the  n^  guess  of  the  value  which  satisfies  the  equation. 


For  C « f[up]  F(uf)  * 0 

CB- 11) 

and 

Fbn) 

up  * uFn neglecting  the  higher  terms  . 

F’  K) 

(B- 12) 
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Using  the  above  formula  we  may  guess  a value  upn  which  satisfies  the 

equation  C = f ^upj  and  calculate  from  Equation  (B-12)  a new  value  of  up 

which  is  a better  fit.  This  new  value  can  then  be  substituted  for  Uc 

rn 

to  improve  the  fit  still  further,  and  so  on.  In  line  720  of  the  program 
the  fit  is  considered  good  enough  if  up  - upn  ^ 10'6. 

Once  the  carrier  concentration  is  known  the  theoretical  trap- 
free  curve  can  be  calculated  from  Equations  (B-l)  to  (B-7),  knowing  the 
oxide  thickness.  This  theoretical  curve  is  then  compared  with  the 
experimental  curve. 

The  curve  fits  are  made  using  the  IMSL  program  package. 

The  program  is  shown  below. 


INTEGER  Ql,  Q2,  P 
REAL  NI,  N2,  LI 

DIMENSION  V (4  0 ) , CAP (40 ) , VX  (4 0 ) , Q (4  0)  ,XK  ( 1 5 ) , Y ( 1 5 ) , 

&C ( 14 , 3 ) 

DIMENSION  WK(630) ,DS (40) ,DSS (1 ) ,US (40) ,CX(40) ,WX(40) 

DIMENSION  UST (40) ,VT(40) ,CAPT (40) ,AB (40) 

READ (11,4)  N2,W0,T,N 

C N2-INITIAL  GUESS  AT  CARRIER  CONC ; WO-OXIDE  THICKNES S ; T-TEMP . 

C N-NO.  OF  DATA  POINTS ; A 1 ,A2 , A3, A4  FIELD  OF  24  CHARACTS  FOR 
C SAMPLE  NAME 

READ (11,4)  (VX(I)  ,1-1 ,N) 

READ (11,4)  (CX(I) ,1-1 ,N) 

READ (11,576)  A1,A2,A3,A4 
576  FORMAT (4A6) 

BB  - 1.0/ ( 1 .005*CX(N)) 

DO  26  I- 1 , N 
AB(I)-CX(I) 

C NORMALISE  CAPACITANCE  VALUES  TO  1.005*MAX  CAP  VALUE 

26  CX(I)  - CX(I )*BB 

4 FORMAT ( ) 

PRINT  11 

11  FORMAT(5X, 'OFFSET  VOLTAGE  IS') 

READ  4, VO 
DO  27  1-1 , N 

27  VX(I)-VX(I)-V0 
J-l 

PRINT  55 

C KX-NO  OF  KNOT  POINTS  FOR  CURVE  FIT 
READ  4 , KX 

55  FORMAT (12X,'TYPE  INTEGER  >2  FOR  NUMBER  OF  KNOTS') 

P-0 

C5-0 

C Vl-Q/KT 

Vl-0.02585*T/300 

C NI-4.9E15(MDE*HDN/M0**2)**.75*EXP(-EG/2KT) 
N1-4.05E14*SQRT(T**3)*EXP(-0.7172/V1) 

PRINT  202 

202  FORMAT (////, 12X, 'A  PROGRAM  THAT  COMPARES  C-V  DATA  TO  THEORY') 
PRINT  274 

274  FORMAT (1  OX, 'AND  COMPUTES  QSS  FROM  THE  VOLTAGE  DIFFERENCES'//) 
PRINT  13,  A1,A2,A3,A4 

13  FORMAT(12X, 'SAMPLE  NUMBER  ',4A6,//) 

CO-KO*EPSO/WOX 

C0-6.28E-9/W0 

C LD-SQRT (KS*EPSO/2Q*NI*KT/Q) 

L 1-1 . 82E3*SQRT (V 1/N 1 ) 

C CD-KS*EPSO/LD 

C1-1.062E-12/L1 
C FERMI  LEVEL-U 1 

U1  - ALOG (N2/N 1 ) 

C CALCULATE  SAMLE  CARRIER  CONC  BY  AN  ITERATIVE  PROCESS  FROM 
C THE  MAXIMIMUM  AND  MINIMUM  CAP  VALUES 


B - 1.062E-12*(I.0/CX(l)  - 1.0)/(C0*L1) 

A - SQRT(8.*Ul-12.+4.*AL0C(Ul-3.)) 

E - B*EXP ( U 1 / 2 ) 

F IS  F (X ) 

F - A + 2.52/ SQRT (U 1 ) - E 
FX  IS  F ' ( X ) 

FX  - (4.*U1-10.)/((U1-3.)*A)-1.26/(U1*SQRT(U1))-0.5*E 
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700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
'6 
810 
820 
830 
840 
850 
860 
870 
' VX  ' 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1))) 
1020 
1030 
2 *F) 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 


C U2-U1-F(X) /F'(X) 

U 2 - Ul  - F/FX 

IF ( AB S (U 1 -U 2 ) - 1 .0E-6*ABS(U1  ) ) 20.20,10 
10  U1  - U 2 

GO  TO  5 

20  N2  - N1*EXP(U1) 

E - Ul*Vl 
C 

PRINT  203,  T,  WO,  CO,  N2,  Nl,  E,  U1 

203  F0RMAT(2X, 'INPUT  AND  CALCULATED  PARAMETERS  ***',//, 

612X, 'TEMPERATURE  -',F6.i,'  DEGREES  K' ,/, 1 2X , 'OXIDE  THICKNESS  - 

4 ,F7.4, ' MICRONS' ,/, 12X, 'OXIDE  CAPACITANCE  -'.1P1E10.3, 

6 ' F/SQ  CM', /,12X, 'DOPING  LEVEL  -'.1P1E9.2,'  PER  CC ' , / , 1 2X , 

&' INTRINS IC  LEVEL  -' , 1 P IE  1 0 . 3 , ' PER  CC ',/, 1 2X ,' FERMI  LEVEL  -' , 
60P1F7.4,'  EV',/,12X, 'NORMALIZED  FERMI  LEVEL  -',F7.3,//> 

PRINT  204,  N 

204  F0RMAT(7X, 'THE  LIST  OF  NORMALIZED  INPUT  C-V  DATA  *',/, 

&1 OX, 'THE  NUMBER  OF  POINTS  IS ' , I 3 , / , l 2X , ' I ' , 7X , 'C  ' , 8X , 'CX' ,8X , 

*,/,  1 IX, '***' ,5X, '****' ,6X, '****' ,6X, '*****' ) 

DO  206  I - 1,  N 

206  PRINT  207,  I,AB(I),  CX(I),  VX(I) 

207  FORMAT (I13,F10.1,F10.3,F10.1) 

C 

C CALCULATE  IDEAL  TRAP  FREE  C-V  CURVE 
C 

C U10  IS  STARTINC  VALUE  OF  US 
U10— .2 

44  DO  300  M-1,31 

IF  (ABS(UIO)  .LT. 0.001 ) U10-0.001 
U-U10/V1 

C F(US,UF)-(2(C0SH(US-UF)-C0SH(UF)+US*SINH(UF) ) )**.5 

F-SQRT ( EXP (U-U 1 )+EXP (U1-U)-EXP(U1)-EXP(-U1 )+U* (EXP(Ul)-EXP(-U 

C CS-CD*SICN(US)*((SINH(US-UF)+SINH(UF))/F(US,UF) ) 

C2-C 1*(U/SQRT(U**2))*(EXP(U-U1)-EXP(U1-U)+EXP(U l)-EXP(-Ul) )/( 

C3-l./(l+(C0/C2)) 

V0-(U+(U/SQRT(U**2))*C1*F/C0)*V1 

C PUT  U AND  V VALUES  IN  ARRAY 

400  UST (M)-U 10 
VT(M)-VO 

IF(U10.EQ.0.001)U 10-0.0 
U10-U 10+0 . 05 

C PUT  CAP  VALUES  IN  ARRAY 
CAPT (M) -C3 
IF (P ) 401,401,300 

C P SETS  LIMIT  ON  THE  MIN  VALUE  OF  CAP.  IT  AVOIDS  POINTS 
C ON  THE  INVERSION  SIDE  OF  THE  CV  CURVE 

401  IF ( CAPT (M) »LE.CX(1 ))  P-M-l 

IF (CAPT (M) .GT .CAPT (M-l ) ) P-M-l 
300  CONTINUE 
253  PRINT  6 

6 F0RMAT(////,2X, 'OUTPUT  DATA  - CUBIC  SPLINE  FIT  ****') 

PRINT  7,  P 

7 FORMAT ( / / , 1 2X , 

4'NUMBER  OF  MATCHING  THEORETICAL  AND  XPMTL  POINTS  IS', 13,//) 

DO  787  01-1, P 
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1260 

1270 

1280 

1290  787 
1300 
1310 
1320 
1 330 
1340 
1350 
1360  85 
1370  C 


Q2-P+1-Q1 
CAP1Q1 )«CAPT (Q  2 ) 
V(Q1)-VT(Q2) 

US !Q1  )«UST!Q2) 

KX1  « KX  - 1 
H - KX  1 

XK(1 ) - 0 . 999  *CX( 1 ) 
XK(KX)  - 1.0 
H - (XK(KX)-XK(l ))/H 
DO  85  1-2, KX1 
XX (X)  - XK(I-l)  + H 


1380  C CURVE  FIT  TO  EXPERIMENTAL  RESULTS 

1390  90  CALL  ICS VKU  (CX , VX  , N , XK  , KX , Y ,C  , 14  , ERROR  , WK , IER  ) 

1400  PRINT  92,  ERROR 

1410  92  FORMAT!//, 12X, 'ERROR  IN  FITTING  THE  DATA  IS ' , IP  1 E8 . 2 , / / ) 
1420  PRINT  94,  (XK(I),I«1,KX) 

1430  94  F0RMAT(12X, 'THE  KNOT  POINTS  ARE ' , / , ( 1 P 5E  1 3 . 3 ) ) 

1440  11  - XK(KX)-XK(KX1  ) 

1450  Y(KX)  - ((C(KX1,3)*H+C(KX1,2))*H+C(KX1,1))*H  +YIKX1) 

1460  C GENERATE  VOLTAGE  VALUES  AT  EXPERIMENTAL  CAP  POINTS 
1470  CALL  ICSEVU(XK,Y,KX,C,14,CAP,WX,P,IER) 

1480  H - CO/1 . 6023E-19 

1490  DO  555  I-l.P 


1500 

L2-P+1 

-I 

1510 

555 

Q!L2) 

- ( V ! I ) - 

WX (I ) )*H 

1520 

PRINT 

1 

1530 

1 

FORMAT 

!6X, 'UST' 

» 

1 2X , 'D V ' 

1540 

DO  2 J 

-1,P 

15  50 

DV-Q ! J 

) / H 

1560 

2 

PRINT 

254 ,UST(J 

) 

,DV 

1570 

254 

FORMAT 

! 1P2E12. 3 

> 

1575 

C CURVE  FIT 

US  VERSUS 

QS 

/) 


1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1655 

1660 

1670 

1680 

1690 

1700 

1710 
1720 
1730 
174  0 
1750 
1760 
1770 
1775 
1780 
1790 
1800 
1810 


H-KX1 

XK(1 )« 1 . 01 *UST  ! I ) 

XK(KX)-1 .05*UST(P) 

H-(XK(KX)-XK(1 ) ) /H 
DO  93  1-2, KX1 
93  XK ( I ) -XK ( I - 1 )+H 

CALL  ICSVKU (UST,Q,P,XK,KX,Y,C, 14, ERROR , WK , I ER ) 

PRINT  92,  ERROR 
C DIFFERENTIATE 

CALL  DCSEVU ( XK , Y , KX , C , 1 4 , US T , DS , P , DS S , P , IE R ) 

PRINT  100 

100  FORMAT!////, 7X, 'OUTPUT  TABLE  OF  QSS  VERSUS  US',//,2X, 

6'I',4X,'US',5X,'INT  QSS ', 4X , 'QSS /EV' , 4X ,' V (THEOR )'  , 4X ,' C ' 
&6X , ' VX' ,6X,'QC',/,1X,'***' ,2X, '****' , 4X, '*******'  , 3X , ' **  * 

63X, '********' , 3X , ' ***' ,4X, '****' ,4X, '****'  ) 

DU  - 0.025 
V(P ) - 0.0 
DO  31  I - 2 , P 
J - P + 1 - I 

31  V(J)  - DU  * !l)S  ( J ) + DS  ! J + l ) ) + V ( J + l ) 

DO  105  I-1,P 

C REFERENCE  ENERGY  VALUES  TO  MID  GAP 
UST (I  ) - UST  ! I ) - E 
L2  - P + 1 - I 
V! I ) - V ! I ) + Q (I > 

CX! I ) - CAP !L 2 ) 


69- 


» » 


1820 

VX(I)  - WX(L2) 

1830 

105 

PRINT  U0,I,UST(I),Q(I) 

1840 

110 

FORMAT ( I 3 , F7 . 3 , IP  2E  1 1 . 3 

184  5 

C CALCULATE  CONSTANT  PART  OF 

1850 

E - 0 

1860 

S - 0.0 

1870 

DO  37  I - 5,P 

1880 

J - P + 4 - I 

1890 

E - E + V(J)*V(J) 

1900 

37 

S - S + V(J) 

1910 

EE  - P - 5 

1920 

SS  - P - 4 

1930 

S - S/SS 

194  0 

E - SQRT((E-SS*S*S)/EE) 

1950 

PRINT  101, S,E 

1960 

101 

F0RMAT(//,12X, 'AVERAGE  i 

1970 

6E11 .3,/) 

1980 

DV  - 1.6023E-19*S/C0 

1990 

PRINT  102,  DV 

2000 

102 

FORMAT (12X,' THE  AVERAGE 

2010 

999 

END 

END  OF  FILE 

QSS 


CONST  QSS 


, 1P1E11.3,'  + /-' 


OFFSET  VOLTAGE  IS  THEN' , F7 . 2 , //) 
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